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ABSTRACT 

G328.4+0.2 is a moderately small (5' x 5') Galactic radio supernova remnant (SNR) at a distance of at 
least 17 kpc that has been long suggested to be Crab-like. Here we report on the detection with ASCA of 
the X-ray emission from the SNR. The X-ray source is faint with an observed flux of (6.0 ± 0.8) x 10~ 13 
erg s _1 cm -2 over the 2-10 keV band. The emission is heavily cut-off at low energies and no flux is 
detected below 2 keV. Spectral analysis confirms that the column density to the source is indeed large, 
Ah ~ 10 23 atoms cm -2 , and consistent with the total column density of hydrogen through the Galaxy 
at this position. Good fits to the spectrum can be obtained for either thermal plasma or nonthermal 
power-law models, although the lack of detected line emission as well as other evidence argues against 
the former interpretation. The power-law index we find, a P = 2.9t°;->, is consistent with other Crab-like 
SNRs. In the radio band G328.4+0.2 is nearly as luminous as the Crab Nebula, yet in the X-ray band 
luminosity it is some 70 times fainter. Nevertheless its inferred soft X-ray band luminosity is greater 
than all but the brightest pulsar-powered synchrotron nebulae and implies that G328.4+0.2 contains a 
rapidly spinning, as yet undetected, pulsar that is losing energy at a rate of ~10 38 erg s . 

Subject headings: ISM: individual (G328.4+0.2, MSH 15-57) 
remnants — X-rays: ISM 



pulsars: general — supernova 



1. INTRODUCTION 

Only about one-quarter of all cataloged Galactic radio- 
emitting supernova remnants (SNRs) (Green 1998) are 
known to emit X-rays based largely on observations by 
the Einstein Observatory and ROSAT. Some, and perhaps 
even many, of the non-detections are likely a result of the 
absorption of the typically soft (i.e., <2 keV) X-ray emis- 
sion from SNRs by the large column density of gas and 
dust in the Galactic plane. However, young remnants, 
like Cassiopeia A, Tycho, and Kepler that show signifi- 
cant thermal emission above 2 keV, and pulsar-powered 
synchrotron nebulae, like the Crab Nebula with hard, fea- 
tureless power-law spectra, should be able to shine through 
the interstellar medium. In fact the entire Galaxy is effec- 
tively transparent to X-rays with energies above about 4 
keV. As part of an effort to identify and study such SNRs, 
we have targeted a sample of small diameter Galactic rem- 
nants with the Advanced Satellite for Cosmology and As- 
trophysics (ASCA), which has high sensitivity and moder- 
ate imaging capabilities in the hard X-ray band. Here we 
report on the discovery of X-ray emission from one of the 
remnants in the sample, G328.4+0.2. 

G328.4+0.2 (MSH 15-57) has been cataloged as a ra- 
dio SNR for quite some time (see, e.g., Clark & Caswell 
1976) and, because of its flat radio spectrum, S v oc i/ -0 - 24 
(Milne 1979), and "filled-center" or plerionic morphology 
(Caswell et al. 1980; Whiteoak & Green 1996), it has long 
been considered as a possible Crab-like remnant. Its dis- 
tance, based on H I absorption studies (Caswell et al. 1975; 
Gaensler, Dickel, & Green 2000 hereafter GDG), is rather 
large, 17 kpc or more, corrected to currently accepted val- 



ues for the rotation curve and distance to the Galactic cen- 
ter. No pulsar has yet been detected (Kaspi et al. 1996) 
and previous searches for X-ray emission have been unsuc- 
cessful (Wilson 1986). 

2. ANALYSIS 

We observed G328.4+0.2 with ASCA (Tanaka, Inouc, 
& Holt 1994) on 17 March 1998. The solid-state imaging 
spectrometer (SIS), which comprises two separate arrays 
of CCD cameras, was operated in 2-CCD mode. The gas 
imaging spectrometer (CIS), consisting of two gas scintilla- 
tion proportional counters, was operated in PH mode with 
the re-assignment of 2 bits of telemetry from pulse-height 
(compressing the spectra from 1024 channels to 256) to 
timing. This, in fact, resulted in no loss of spectral res- 
olution since the GIS spectra are re-binned to this reso- 
lution before spectral analysis anyway. Rise-time infor- 
mation was used for suppressing background events. The 
data were screened using standard criteria that resulted 
in a combined exposure of 32900 s for the SIS detectors 
and 38950 s for the GIS detectors. The data from each de- 
tector were separately reduced and like observations were 
combined. 

Images (see Figure 1) were made in two energy bands, 
which were chosen to optimize the signal-to-noise ratio 
of the detected sources. In each of the four sets of data 
(i.e., the two separate GIS and SIS detectors) two sources 
were apparent. One, a soft unresolved object, was de- 
tected with a statistical significance of 8. Oct in the SIS and 
7.5ct in the GIS. This source is coincident with a previ- 
ously known X-ray source from the Einstein Observatory, 
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Fig. 1. — (a) X-ray image of the SNR G328.4+0.2 and environs from a 32900 s observation with the ASCA SIS in soft and hard X-ray 
bands as indicated. Plus signs denote the location of the radio SNR (the source to the west) and a previously known soft X-ray source. The 
data were smoothed with a 1' (a) gaussian. In the soft band image, the greyscale runs linearly from 1.87 X 10 — 4 SIS counts s _1 arcmin -2 
(3 a above background) to 5.76 X 10 — 4 SIS counts s — 1 arcmin - 2 (the peak emission). In the hard band image, the greyscale runs linearly 



from 3.00 X 10" 4 SIS counts s" 



(3 a above background) to 4.55 X 10 4 SIS counts s 



(the peak emission). In each panel 



contours are plotted at linear intervals from the minimum value to 90% of the peak. The dotted contour indicates the edge of the detector. 
The approximate size of the point-sprcad-function is indicated in the lower right. 




Fig. 1. — (b) Same as (a) except for the following. The data are from a 38950 s observation with the ASCA GIS. In the soft band image, 
the greyscale runs linearly from 1.19 X 10 — 4 GIS counts s _1 arcmin" 2 to 2.92 X 10 — 4 GIS counts s — 1 arcmin -2 . In the hard band image, the 
greyscale runs linearly from 2.57 X 10 -4 GIS counts s _1 arcmin -2 to 5.01 X 10~ 4 SIS counts s _1 arcmin -2 . 



2E 1552.8-5316, (Wilson 1986; Harris et al. 1990) and the 
ROSAT all sky survey, 1RXS J155644.7-532441 (Voges et 
al. 1999). We will not discuss this source any further, other 
than to note that it appears only in the band below 2 keV 
and is therefore likely to be a relatively nearby Galactic 
object. 

In the band above 2 keV, at a position consistent with 



that of G328.4+0.2, we detect a new X-ray source (des- 
ignated AX J1555.5— 5318) with a statistical significance 
of 8.4cr in the SIS and 8.5cr in the GIS. There are roughly 
167 ± 20 SIS events and 222 ± 26 GIS events from the new 
source, which is unresolved at the spatial resolution af- 
forded by the ASCA telescopes (i.e., roughly 3' half power 
diameter). There are no significant X-ray sources else- 
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Fig. 2. — X-ray spectra of G328.4+0.2 from the ASCA SIS and GIS data as indicated. The best-fit power-law spectral model is shown as 
the histogram. 



where in the fields of view of either the SIS or GIS other 
than the two we introduce above. 

For completeness we carried out a search for coherent 
pulsations in AX J1555. 5-5318 using the 2-10 keV GIS 
data. As a consequence of the re-assignment of some 
telemetry to timing information, as mentioned above, the 
time resolution of the GIS data was 0.015625 s and 0.125 
s during the high and medium bit rate portions of the 
observation. These data were separately Fourier trans- 
formed (after correcting photon arrival times to the solar 
system barycenter) and the resulting power spectra were 
examined for significant peaks. The periods searched were 
0.03125 s to 1024 s for the high time resolution data and 
0.25 s to 4096 s for the medium time resolution data. No 
significant pulsed signal was detected, which is not sur- 
prising given the small number of source photons. We set 
a broad limit on the pulsed fraction of < 50%. 

Source and background spectra were extracted from 
both the SIS and GIS. The source regions were circles cen- 
tered on the new X-ray source position with radii, chosen 
to optimize the signal-to-noise ratio of the extracted spec- 
tra, of 3!2 (SIS) and 3'.9 (GIS). The background region for 
the SIS was taken from the apparently blank portion of 
the field of view surrounding the source toward the NW, 
SW, and SE (see Fig 1). For the GIS background we used 
an annular region centered on the position of the optical 
axis with inner and outer radii chosen so that the annulus 
just contained the circular source region. This was done 
to ensure that the background region contained emission 
that had experienced the same off-axis telescope properties 
(e.g., vignetting) as the source emission. Of course, both 
the source region and a region encompassing the X-ray 
emission from the other soft X-ray source were excluded 
from the GIS background annulus. 

Local background regions were used because of concerns 
that diffuse emission from the well-known Galactic ridge 
emission (Worrall et al. 1982) might be present in the field 



of view. In order to address this point the GIS spectra of 
the background regions were themselves studied using the 
new (point-source removed) blank-field background files 
from high Galactic latitudes. There is a clear excess of X- 
ray emission in the background regions near G328.4+0.2 
that can be fit well by an absorbed thermal plasma model 
(Raymond & Smith 1977) with kT ~ 6 - 7 keV and a 
surface brightness of 1.3 x 10~ 4 counts s _1 arcmin~ 2 or 
cquivalcntly (in terms of flux over the 2 10 keV band) 
1.1 x 10~ 14 ergs s _1 cm~ 2 arcmin~ 2 . This is compara- 
ble to the brightness level expected for the ridge near this 
position. Our use of a local background is, therefore, not 
only justified but essential. 

Response files for spectral fitting were generated em- 
ploying standard software tools in the usual manner for 
ASCA. The combined SIS and GIS spectral data were re- 
grouped so that each spectral bin contained a minimum 
of 25 events in order to allow use of x 2 as the figure- 
of-merit function for spectral fitting. Given the limited 
statistical quality of the data, only two simple spectral 
models were explored: an absorbed power-law and an ab- 
sorbed thermal plasma model (Raymond & Smith 1977). 
The spectral data and power-law model are shown in Fig- 
ure 2 and numerical values for the best-fit parameters are 
given in Table 1 . A slightly better fit was obtained for the 
power-law model, although both models were statistically 
acceptable. The observed 2-10 keV band flux of the source 
is (6.0 ± 0.8) x 10~ 13 erg s" 1 cm" 2 . 

3. DISCUSSION 

The large absorption we derive from both sets of spectral 
fits is fully consistent with the total Galactic column of hy- 
drogen along the line-of-sight to G328.4+0.2, A H ~ 10 23 
atoms cm -2 . The column can be decomposed into three 
parts: (1) a neutral hydrogen column of Am ~3x 10 22 
atoms cm~ 2 (Wilson 1986), (2) an ionized hydrogen col- 
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Table 1 

Spectral Model Fits for G328.4+0.2 



Parameter Value and Error (1 a) 
Power-law Model 

JV H (H atoms cm" 2 ) 9.1±|"2 x 10 22 

O Q+0-9 

"P 2 - 9 -0.8 

F E (lkeV) (photon s" 1 cm" 2 kcV" 1 ) 1.8±f"^ x 10~ 3 

X 2 /d.o.f 37^6/48 

Thermal Plasma Model 

N H (H atoms cur 2 ) 7.4± 2 ' 2 x 10 22 

kT (keV) 4.b±f"| 

Fractional abundance (©) <0.68 (90% C.L.) 

Emission measure (cm" 5 ) 1.4+J^ x 10 11 

X 2 /d.o.f 38'.2/47 



umn of JVhh ~ 3 x 10 22 atoms cm -2 estimated from the 
free-electron model of Taylor and Cordes (1993), and (3) 
a molecular hydrogen column of Nn 2 ~ 2 x 10 22 molecules 
cm -2 based on the CO data from Bronfman et al. 1989 and 
Bitran et al. 1997 (as reported by Slane et al. 1999). The 
agreement between these measures of the total hydrogen 
column supports the large distance of 17 kpc to this rem- 
nant derived from the H I absorption studies mentioned 
above. 

3.1. Thermal Interpretation 

The spectral fits themselves do not preclude a thermal 
interpretation for the X-ray emission from G328.4+0.2. 
However there are two main difficulties with the interpre- 
tation that cause us to disfavor it. First and foremost is 
the combination of a moderately high temperature and 
a lack of significant Fe-line emission. Mean plasma tem- 
peratures of greater than 2 keV tend to be the exception 
rather than the rule for SNRs and in all cases to our knowl- 
edge are accompanied by Fe Ka line emission at ~6.7 keV. 
The second problem arises from the inferred evolutionary 
state of the remnant. It is possible to determine the age 
and explosion energy from the remnant's measured size, 
temperature, and emission measure in the context of sim- 
ilarity solutions for SNRs in the adiabatic phase of evolu- 
tion (see, for example, Hughes, Hayashi, & Koyama 1998). 
The results are sensitive to assumptions about the state of 
equilibration between electrons and ions. If one assumes 
that the temperatures of the two species are fully equili- 
brated everywhere throughout the remnant, then the age 
and explosion energy of G328.4+0.2 are estimated to be 
2000-4000 yr and (0.9 - 1.8) x 10 51 erg, using values and 
uncertainty ranges from Table 1 and a total size of 2!5 in 
radius (GDG). Although these results are nicely consistent 
with expected values, various lines of evidence now suggest 
that full equilibration of electron and ion temperatures at 
supernova shock fronts is unlikely to be the correct sce- 
nario (Laming et al. 1996, Hughes et al. 1998). If one 
relaxes the assumption of full equilibration in favor of a 
scenario where Coloumb collisions regulate the exchange 
of energy between the species, then there are no plausible 



solutions for the age and explosion energy of G328.4+0.2 
(i.e., the inferred explosion energies are at least two or- 
ders of magnitude larger than the canonical value of 10 51 
erg). In summary, we consider it extremely unlikely that 
G328.4+0.2 is a young, thermally-emitting SNR. 

3.2. Nonthermal Interpretation 

Based on the radio properties of this SNR, specifically 
the filled-center morphology and high degree of linear po- 
larization of the nebula (^20%; GDG), a nonthermal in- 
terpretation of the X-ray emission is most probable. And 
indeed the power-law index we find is consistent with that 
of other plerionic SNRs, of which the brightest and best 
known is the Crab Nebula with ap = 2.1 (Toor & Seward 
1974). The nebular luminosity of such sources is powered 
by the spin-down energy loss of a rapidly-rotating young 
neutron star and the emission mechanism is synchrotron 
radiation from relativistic electrons in the presence of the 
magnetic field of the nebula. In some cases, like the Crab, 
pulsed radiation from the central compact star is observed, 
although in other cases it is not, perhaps due to unfavor- 
able viewing geometry of the presumably beamed radia- 
tion from the pulsar. However, emission from the syn- 
chrotron nebula is isotropic and therefore, when detected, 
a synchrotron nebula is believed to provide clear evidence 
for the presence of a young pulsar. 

The unabsorbed soft X-ray luminosity (0.2-4 keV band) 
of G328.4+0.2islog(L x /ergs- 1 ) = 35.6±J;?. This value is 
larger than all known pulsar-powered synchrotron nebula 
with the exception of the Crab (Seward & Wang 1988), and 
two SNRs in the Large Magellanic Cloud: E0540-69.3 (Se- 
ward & Wang 1988) and N157B (Wang & Gotthelf 1998). 
It is comparable to the luminosities of the Crab-like cores 
of G29.7-0.3 (Blanton & Helfand 1996), G21.5-0.9 (Slane 
et al. 2000), and MSH 15-52 (Seward & Wang 1988). 
From the empirical relationship between the soft X-ray lu- 
minosity of the synchrotron nebula, Lx and the pulsar's 
rate of loss of rotational energy, E (Seward & Wang 1988), 
we determine that the (still undiscovered) pulsar that must 
be powering the synchrotron emission is losing rotational 
energy at the rate of E ~ 10 37 - 2 x 10 38 ergs _1 . GDG use 
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the radio luminosity to estimate a rotational energy loss 
rate of several 10 38 ergs _1 . Given the large scatter in the 
correlations between E and luminosity in each band, the 
differences between these estimates are unlikely to be sig- 
nificant. Indeed, focusing on the more interesting broad 
agreement between them, we conclude that G328.4+0.2 
harbors a rapidly spinning pulsar that is losing rotational 
energy at a rate of ~10 38 ergs _1 . 

GDG argue that the size and radio luminosity of 
G328.4+0.2 most closely resemble those of N157B, the 
SNR that contains the most rapidly spinning known pul- 
sar with a period of 0.016 s (PSR 0537-6910; Marshall 
et al. 1998). This is also the case for the X-ray luminos- 
ity and for the ratio of X-ray to radio luminosity. The 
X-ray luminosity of the power-law component in N157B is 
log (Lx /ergs -1 ) = 36.7, at the upper end of our allowed 
range for G328.4+0.2. Although the ratio of X-ray to ra- 
dio luminosity is low for G328.4+0.2, L X /L R = 0.3 - 25, 
nevertheless it is an excellent match to that of N157B 
(Lx I Lr ~ 10) . In contrast Lx / Lr arc 130 and 800 for the 
Crab and E0540— 69.3, respectively. According to Wang 
& Gotthelf (1998) the X-ray synchrotron nebula in N157B 
is some 5 x 7 pc in size and is surrounded by a fainter 
diffuse component of thermal X-ray emission from the SN 
blast wave. Because of the broad point-spread-function 
of ASCA, we are unable to place interesting limits on the 
size of the X-ray synchrotron nebula in G328.4+0.2. This 
must await more sensitive and higher spatial resolution 
X-ray observations. Unfortunately, the large absorption 
toward G328.4+0.2 may forever preclude the detection of 
soft thermal X-ray emission (with kT ~ 1 keV) from any 
outer shell or blast wave in G328.4+0.2. We note that 
there is no evidence for a shell from the high spatial reso- 



lution radio image (GDG). 

The broadband spectrum of G328.4+0.2 from radio to 
X-ray frequencies is shown in Figure 3. The point at which 
the extrapolated best-fit radio and X-ray power-laws inter- 
sect, the so-called break frequency, is v-q ~ 3 x 10 14 Hz. 
The large uncertainty on the measured X-ray spectral in- 
dex results in a huge range for v-q (10 11 Hz - 6 x 10 15 Hz) 
that provides little more than consistency with values from 
other SNRs, e.g., the Crab with vr ~ 10 13 Hz. Compared 
to the Crab, G328.4+0.2 is nearly as luminous in the ra- 
dio band (GDG), while it is some 70 times less luminous 
in X-rays. 

Synchrotron radiation theory (e.g., Ginzburg & Sy- 
rovatskii 1965) provides a means for estimating from the 
observed spectrum the amount of energy in relativistic par- 
ticles and magnetic field in the nebula. One needs an ad- 
ditional constraint which we take to be the condition that 
the total energy is at a minimum. In this case the mag- 
netic field energy is 3/4 of the energy in particles. We 
further assume that the energy in all relativistic particles 
is twice the energy in the electron component alone. We 
integrate the broadband spectrum (Fig. 3) from low fre- 
quency to i>b using the radio spectral index and from vb 
to 2.4 x 10 18 Hz (10 keV) using the X-ray spectral index. 
For our derived quantities we show explicitly the depen- 
dencies on distance, D, angular radius of the nebula, 9, the 
ratio of the energy in magnetic field to that in particles, 
K m , and the ratio of energy in all particles to electrons, 
k t . We use a value, 9 = 1', for the size of the nebula 
that corresponds to the bright central portion (GDG). 
The value we derive for the nebular magnetic field is 
B = 365( Km /O.75) 2 / 7 («; r /2) 2 / 7 (Z?/17kpc)- 2 / 7 (0/T)- 6/7 
[iG and the energy in relativistic electrons is W e = 5 x 
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Fig. 3. — Broadband spectrum of G328.4+0.2 from radio to X-ray frequencies. Radio data points come from Caswell et al. (1975) (1.42 
GHz); Clark & Caswell (1976) (0.408 GHz and 5 GHz). Tateyama, Sabalisck, & Abraham (1988) (5 GHz and 22 GHz); and Whiteoak & 
Green (1996) (0.843 GHz). The dashed curves show the extrapolation of the best-fit radio and X-ray power-law models to the point where 
they intersect at the so-called break frequency, which is ~3 X 10 14 Hz (albeit with large errors). 
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10 49 (L»/17kpc) 2 ( J B/365^G)- 3/2 erg, with roughly equal 
amounts of energy in the populations of electrons above 
and below the break frequency. Our estimates of the 
above quantities are reduced to B ~ 220 and W e ~ 
2 x 10 49 erg if we use the lower limit on the break frequency 
(10 11 Hz) that comes from the 1-a lower limit on the X- 
ray spectral index. When the total energy in the nebula 
is computed, viz., W = W cr + Wb — «v(l + n m )W e , even 
the lower limit value is still rather high, W = 7 x 10 49 erg. 

Although the energy we see currently in the synchrotron 
nebula must have come from the initial rotational energy 
of the pulsar, Eq = 2tt 2 IP^ 2 , it is also clear that W can 
represent only a fraction of Eq. Some of the initial ro- 
tational energy of the pulsar has gone into radiation and 
the expansion of the nebula and some still remains in the 
form of the current spin of the pulsar. Calculating the 
precise partition of Eq into the various forms that it may 
take is a challenging task and, given the limited quality of 
the current data, is beyond the scope of this work. Nev- 
ertheless by equating W to Eq we can derive a strong 
upper limit to the initial spin period of the pulsar. Us- 
ing I = 10 45 g cm 2 for the moment of inertia, the limit we 
obtain is Pq < 0.017 s. This limit is fully consistent with 
the current spin period of PSR 0537—6910 and adds fur- 
ther weight to the preceding arguments, as well as those 
in GDG, that the closest known analogue to G328.4+0.2 
is the Crab-like SNR N157B. 

4. CONCLUSIONS 

Using ASCA we have discovered hard X-ray emission (2- 
10 keV band) from the Galactic radio SNR G328.4+0.2. 



The new X-ray spectral data can be equally well fit by 
highly absorbed (Nr ~ 10 23 atoms cm -2 ) thermal plasma 
or power-law models. When the inferred astrophysical 
properties of the remnant under each scenario are exam- 
ined, along with additional constraints from recent radio 
observations, we conclude that G328.4+0.2 is a Crab-like 
SNR dominated by synchrotron emission in the radio and 
X-ray bands. It is powered by the spin-down energy of a 
central, yet undetected, pulsar that is currently losing en- 
ergy at a rate of ~10 38 erg s _1 . In order to account for the 
total energy in particles and magnetic field observed in the 
synchrotron nebula at the present time, we find that the 
initial spin period of the pulsar must have been no more 
than 0.017 s. In many respects G328.4+0.2 bears a partic- 
ularly close resemblance to the Crab-like SNR N157B in 
the Large Magellanic Cloud, the remnant containing the 
most rapidly rotating pulsar known, PSR 0537—6910. The 
discovery of its hidden pulsar is clearly an important next 
step in furthering our understanding of G328.4+0.2, yet 
the faintness of the remnant in the X-ray band (6 x 10 -13 
erg s _1 cm -2 over 2-10 keV) will make this search a chal- 
lenging task. 

This work made use of the SIMBAD Astronomical 
Database maintained by the Centre de Donnees As- 
tronomiques de Strasbourg. Useful conversations with 
Cara Rakowski and Brian Gaensler are greatly acknowl- 
edged. We appreciate Monique Arnaud's support and hos- 
pitality during the course of this project. Partial support 
was provided by NASA grant NAG5-6420 (JPH) and con- 
tract NAS8-39073 (POS, PPP). 
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